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DESCRIPTION 



Curve Generating Apparatus and Curve Generating Method, and Program 
Recording Medium 

Technical Field 

This invention relates to a curve generating apparatus and a curve generating 
method used in CAD (Computer Aided Design), CG (Computer Graphics), etc. which 
edit curve indicating contour shape of object included within picture image on the 
computer, and a program recording medium adapted so that program relating to the 
curve generating method is recorded. 

Background Art 

Hitherto, there are many instances where contour extraction processing for 
extracting contour shape of object included in plural picture image data constituting 
moving picture image are used in the field of picture image processing, e.g., CAD, CG 
technologies, etc. Among these picture image processing, when picture image 
synthesis processing is carried out, there is the processing which generates key signals 
from contour shape of the object, and it is important to generate precise or accurate 
key signals, i.e., precise contour shape. In this case, the key signal is information for 
cutting off the area of the foreground object to be synthesized and is also called mask. 
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The contour extraction processing is required to have interactive operability as in the 
case where operator gives in detail information such as contour position or direction, 
etc. with respect to object. As such contour extraction processing, there are processing 
to designate plural points on the contour of object to prepare shapes between 
respective points, processing to operate control points of parametric curve indicating 
counter shape, and processing to directly input contour shape by mouse, etc. As 
examples of these contour extraction processing, there are known "Intelligent Scissors 
for Image Composition" (Eric N. Mortensen and William A. Barrett, Computer 
Graphics Proceedings, Annual Conference Series, 1995, ACM SIGGRAPH, pp. 191- 
198), "Image contour detection method" (Japanese Patent Application Laid Open No. 
152481/1992 publication), and "Cutting mask preparation method and cutting mask 
preparation apparatus" (Japanese Patent Application Laid Open No. 254854/1992 
publication). 

Moreover, in the contour extraction processing, it is required to extract precise 
contour shape with respect to all picture images constituting moving picture. For 
example, in cinemas or television picture images, several hundreds of key signals are 
required for generating images of several seconds, and the quantity of processing is 
vast. For this reason, the contour extraction processing used in the field of cinema, 
etc. is desired (required) that more precise contour shape can be prepared by lesser 
number of operations. This results from the fact that, in the conventional contour 
extraction processing, processing for preparing contour shapes with respect to 
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respective frames constituting moving picture is the processing which requires 
extremely much time and labor. 

In view of the above, as a method of preparing contour shape by more simple 
procedure, there was proposed a method in which, on the premise that, in moving 
picture images having continuous property, contour shape at the start frame changes 
with passage of time so that it changes to contour shape at the end frame, intermediate 
shape of contour at the intermediate frame is prepared by interpolation processing 
from contour shape at the start frame and contour shape at the end frame. 

In accordance with such an interpolation processing, as shown in FIG. 1 , plural 
points constituting contour shape are caused to have one-to-one correspondence 
relationship at contour shape of the start frame and contour shape of the end frame. 
Namely, when it is assumed that curve indicating contour shape at the start frame (time 
t =0) is curve C 0 , curve indicating contour shape at the end frame (time t=l) is curve 
C 1? and curve changes in point of time, it can be considered that curve C 0 at time t=0 
changes to curve Q at time t=l. In view of the above, when time t at which 
interpolation of two curves is assumed to be T, and coordinate at the point on the curve 
C 0 is put by A and coordinate at the point on the curve Q corresponding to the 
coordinate A is put by B, coordinate C at curve C x which takes required intermediate 
form can be represented as follows. 

C= T- A + (l-T)-B 
In the conventional interpolation processing, contour curve C x at the intermediate 
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frame (t=T) was prepared from train of points determined by such interpolating 
calculation. 

Further, it is assumed that contour shapes with respect to respective frames are 
constituted by plural third-order Bezier curve. The third-order Bezier curve is defined 
by 

Q( t )=((l_t) 3 )M+3t«^ t 1) (0<t* 1) 

as indicated by the article "Computer Graphics PRINCIPLE AND PRACTICE 
SECOND EDITION in C" (Foley, van Dam, Feiner, Hughes, ADDISON WESLEY, 
1996, ISBN 0-201-84840-6). In this formula, M, N, O, P are coordinates indicating 
points on the two-dimensional plane wherein M and P are end points of line segment, 
and N, O are control points. An example of line segment of the third order Bezier 
curve represented by the above-mentioned formula is shown in FIG. 2. In FIG. 2, 
curve indicated by solid line is the third-order Bezier curve, wherein M, P indicated 
by points of black circle are end points, and points N, O indicated by white circle are 
control points. 

In accordance with the formula which defines the above-described Bezier curve, 
locus of the third-order Bezier curve, i.e., shape is represented by coordinates of end 
points and control points indicated by points M, N, O, P and change of time t in FIG. 

2. 

When constituting unit of the third-order Bezier curve shown in FIG. 2 is 
assumed to be segment, it can be said that single curve indicating contour shape is 
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ordinarily composed of plural segments as shown in FIG. 3. In this case, end points 
of adjacent segments are shared so that one curve is constituted as the entirety. 

In the above-described conventional contour extraction processing, in the case 
where intermediate shape is determined from contour shapes of frames before and 
after in point of time, interpolation was carried out in the state where control points R 
and end points S constituting segments in FIG. 3 are caused to have one-to-one 
correspondence relationship. 

However, at frames before and after in point of time, there are instances where 
contour shape does not linearly move. 

In the above-mentioned method, in the case where contour shape move non- 
linearly, shape different from contour shape was prepared. 

In order to avoid such problems, it is conceivable to carry out work to divide 
interval in point of time between the time of start and the end time of end into 
narrower ranges to carry out interpolation processing of shape for a second time, or 
to manually modify the shape, etc. For this reason, there took place the problem that 
labour of user is required in preparing shape, and there takes place fine unevenness of 
shape because shape discontinuously changing is prepared in the time axis direction 
owing to manual modification. 

Disclosure of the Invention 

This invention has been made in view of actual circumstances as described 
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above, and its object is to provide a curve generating apparatus and a curve generating 
method in which, in the processing for preparing intermediate shape of two curves, 
even in the case where contour shape does not linearly move, it is possible to precisely 
pursue correspondence points on the contour, thus to prepare contour shapes with 
higher precision. 

For this reason, the curve generating apparatus according to this invention is 
directed to a curve generating apparatus adapted for generating, from curve to be 
generated of first frame, curve at second frame, the apparatus comprising: 
correspondence point detecting means for determining, by using curve to be generated 
of the first frame, correspondence points at the second frame corresponding to 
reference correspondence points set on the curve to be generated of the first frame; 
and curve generating means for generating curve of the second frame passing through 
correspondence points detected by the correspondence point detecting means.. 

Moreover, the curve generating method according to this invention is directed 
to a curve generating method for generating, from curve to be generated of the first 
frame, curve at second frame, the method comprising: a correspondence point 
detection step of determining, by using the curve to be generated of the first frame, 
correspondence points at the second frame corresponding to reference correspondence 
points set on the curve to be generated of the first frame; and a curve generation step 
of generating curve of the second frame passing through the correspondence points 
detected by the correspondence detection step. 
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In addition, the program recording medium according to this invention is 
directed to a program providing medium for providing program relating to curve 
generation processing for generating, from curve to be generated of first frame, curve 
at second frame, wherein there is recorded program comprising: a correspondence 
point detecting step of deteniiining, by using the curve to be generated of the first 
frame, correspondence points at the second frame corresponding to reference 
correspondence points set on the curve to be generated of the first frame; and a curve 
generation step of generating curve of the second frame passing through the 
correspondence points detected by the correspondence point detection step. 

Brief Description of the Drawings 

FIG. 1 is a view for explaining a method of preparing shape of intermediate 
frame by interpolation processing from start frame and end frame. 

FIG. 2 is a view showing an example of line segment of third-order Beizer 

curve. 

FIG. 3 is a view showing single curve indicating contour shape as plural 
segments. 

FIG. 4 is a block diagram showing configuration of a curve generating apparatus 
to which this invention is applied. 

FIG. 5 is a view for explaining moving picture consisting of plural frames to be 
processed at the above-mentioned curve generating apparatus. 




8 



FIG. 6 is a flowchart showing processing procedure in software program 
relating to a curve generating method of this invention. 

FIG. 7 is a view for explaining a method of carrying out pursuit by using 
absolute value sum of differences of picture images as pursuing method for picture 
image. 

FIG. 8 is a view for explaining first contour curve A and second contour curve 
B used when intermediate curve is generated by the curve generating apparatus. 

FIG. 9 is a view showing intermediate curve generated by the curve generating 
apparatus from first contour curve A and second contour curve B. 

FIG. 10 is a flowchart showing processing when intermediate curve is generated 
from the first contour curve A and the second contour curve B. 

FIG. 1 1 is a view for explaining that sampling points are set to determine 
sampling interval. 

FIG. 12 is a flowchart showing processing when sampling intervals with respect 
to first contour curve A and second contour curve B are determined. 

FIG. 13 is a view showing first contour curve A and second contour curve B 
from which sampling interval is obtained. 

FIG. 14 is a flowchart showing processing when resampling processing is 
carried out. 

FIG. 15 is a flowchart showing processing when intermediate curve indicating 
intermediate shape is generated. 
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FIG. 16 is a flowchart for explaining processing procedure for deforming 
interpolated shape in correspondence with reference correspondence points. 

FIG. 17 is a view showing a more practical example of processing up to step 

542 of the flowchart shown in FIG. 16. 

FIG. 18 is a view showing a more practical example of processing at the step 

543 and steps subsequent thereto of the flowchart shown in FIG. 16. 

FIG. 19 is a flowchart for explaining operation of another more practical 
example of this invention. 

FIG. 20 is a block diagram showing an example of the configuration of the main 
part of curve generating means for generating curve by the more practical example of 
path search method. 

FIG. 21 is a flowchart showing basic processing procedure of the path search 
method. 

FIG. 22 is a view for explaining processing for carrying out path search on the 
basis of the path search method. 

FIG. 23 is a view for explaining internal configuration of path search processing 
optimization means. 

FIG. 24 is a flowchart for explaining processing procedure in the path search 
processing optimization means. 

FIG. 25 is a view for explaining processing relating to passing cost calculation 
parameter of the path search processing optimization means. 
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FIG. 26 is a flowchart for explaining processing procedure in the path search 
means. 

FIG. 27 is a view for explaining, in detail, processing relating to passing cost 
calculation parameter utilizing selective gradient detection. 

FIG. 28 is a flowchart for explaining processing procedure in passing cost 
optimization means. 

FIG. 29 A is a view showing one example of passing cost map of path search and 
FIG. 29B is a view showing process where minimum cost path calculation means 
determines path when passing cost map is obtained. 

FIG. 30 is a view showing one example of the shortest path search algorithm 
using dynamic programming. 

FIG. 3 1 is a view for explaining in detail processing relating to calculation range 
optimization of path search processing optimization means. 

FIG. 32 is a flowchart for explaining processing procedure in calculation range 
optimization means. 

FIG. 33 is a block diagram showing an example of the configuration of point 
train generating means. 

FIG. 34 is a flowchart for explaining processing procedure in the point train 

generating means. 

FIG. 35 is a block diagram showing the configuration of a more practical 
example of optimum color projection axis calculating unit for determining color 
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projection axis. 

FIG. 36 is a view for explaining the operation principle of a more practical 
example of the optimum color projection axis calculating unit. 

FIGS. 37 A and 37B are views for explaining the operation principle of another 
more practical example of the optimum color projection axis calculating unit. 

Best Mode for Carrying Out the Invention 

Preferred embodiments of this invention will now be described in detail with 
reference to the attached drawings. 

A curve generating apparatus 1 to which this inventionb is applied is caused to 
be of configuration as shown in FIG. 4, for example. This curve generating apparatus 
1 is adapted so that a program memory 3, a data memory 4, a frame memory 5, an 
operation input unit 6, and an external storage (memory) unit 7 are connected to an 
arithmetic processing unit 2 through a bus 8. 

In the program memory 3, there is recorded program for executing a curve 
generating method according to this invention. Namely, the program memory 3 
records therein, as software program, the curve generating method comprising a 
correspondence point detection step of determining, by using curve to be generated or 
processed of the first frame, correspondence points at second frame corresponding to 
reference correspondence points set on the curve to be generated of the first frame, 
and a curve generation step of generating curve of the second frame passing through 
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correspondence points detected by the correspondence point detection step. 

The operation input unit 6 is comprised of, e.g., keyboard and/or mouse, etc. 
and serves to generate an operation input signal as the result of the fact that it is 
operated by user. This operation input unit 6 outputs the operation input signal to the 
arithmetic processing unit 2. 

The arithmetic processing unit 2 generates control signals with respect to the 
above-described respective units in accordance with operation input signals from the 
operation input unit 6. Further, the arithmetic processing unit 2 outputs control signals 
thus generated to respective units to thereby cany out processing corresponding to the 
operation input signals. Namely, the arithmetic processing unit 2 conducts such a 
control to read out in sequence the software programs from the program memory 3 to 
execute them to generate contour curve indicating contours of curves to be generated 
included in respective frames by using picture data indicating moving picture stored 
in the external recording unit 7. 

Picture image data processed by the software* program from the program 
memory 3 is stored in the external storage unit 7 as described above. The external 
storage unit 7 comprises recording medium, e.g., optical disc, etc. and outputs, in 
accordance with control signal from the arithmetic processing unit 2, picture data 
indicating moving picture recorded on the recording medium to the arithmetic 
processing unit 2, the data memory 4 and the frame memory 5. In addition, the 
external storage (memory) unit 7 records, onto the recording medium, other picture 
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data delivered from the external. 

As shown in FIG. 5, for example, picture image data stored in the external 
storage unit 7 consists of plural frames Frl, Fr2, Fr3, Fr4, Fr5, Fr6 and Fr7 to 
constitute moving picture. In this example, frame Frl is start frame and end frame 
with respect to the frame Frl is frame Fr3. Namely, this curve generating apparatus 
1 executes the program relating to the curve generating method of this invention to 
thereby generate intermediate curve with respect to frame Fr2 which is intermediate 
frame by using first contour curve of data to be generated included in the frame Frl 
and second contour curve of data to be generated included in the frame Fr3. 

In the data memory 4 of FIG. 4, there is stored contour curve indicating contour 
of intermediate frame generated by the above-described arithmetic processing unit 2 
for executing the program, etc. In addition, into this data memory 4, data in the middle 
of arithmetic processing by the arithmetic processing unit 2 is also stored. 

Picture image data corresponding to one frame from the data memory 4 or the 
external storage (memory) unit 7, etc. is inputted to the frame memory 5 in accordance 
with control signal of the arithmetic processing unit 2. Picture image data which has 
been read out by control signal from, e.g., arithmetic processing unit 2 or contour 
curve indicating intermediate frame generated by the program, etc. is inputted to the 
frame memory 5. Further, this frame memory 5 is connected to a display unit 9 for 
displaying picture image and/or processing content, etc. to user, and outputs picture 
image data, etc. inputted in accordance with control signal to the display unit 9 to 
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thereby allow the display unit 9 to display thereon picture image that the above- 
mentioned picture image data indicates. 

The program to which the curve generating method of this invention is applied 
is operated by the processing procedure shown in FIG. 6. 

At step S 1 0 1 , user gives, as input, two curves at the time of start and at the time 
of end and reference correspondence points with respect to respective curves. In this 
case, two curves are respectively referred to as curves A and B. It is required that the 
numbers of reference correspondence points are the same with respect to the curves 
A and B. Then, at step S 102, pursuit operation of reference correspondence points is 
carried out with respect to picture images successive in the time axis direction. As a 
result, positions of reference correspondence points at respective intermediate stages 
can be determined. The detail thereof will be described later. 

Then, at step S103, interpolation of shape is carried out. In order to cany out 
shape interpolation, such an approach is employed to carry out interpolation 
calculation of coordinates of point trains on two curves that user designates at the step 
S101 to thereby determine coordinates of points of intermediate shapes. Although 
there maybe employed any transformation (conversion) to provide intermediate shape 
of two shapes, there is employed in this example the technology disclosed in the 
Japanese Patent Application No. 251492/1998 by the applicant of this application. 
The detail thereof will be described later. 

Then, at step SI 04, the interpolated shape determined at the step SI 03 is 
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deformed in correspondence with reference correspondence points determined by 
pursuit operation at the step SI 02. 

Accordingly, by executing the program, it can be equivalently said that the 
arithmetic processing unit 2 comprises a correspondence point detecting section for 
determining, by using curve to be generated or processed of the first frame, 
correspondence points at intermediate frame corresponding to reference 
correspondence point set on the curve to be generated of frame at the time of start, and 
a curve generating section for generating curve of intermediate frame passing through 
correspondence point detected by the correspondence point detecting section. 

Particularly, the above-mentioned correspondence point detecting section 
carries out picture image pursuit operation at the step SI 02 to detect where reference 
correspondence point on curve to be generated or processed of frame at the time of 
start is positioned at the intermediate frame. 

In addition, the curve generating section generates curve passing through 
correspondence points detected by the correspondence point detecting section through 
the steps S103 and S104. In more practical sense, when frame at the time of end is 
assumed to be the third frame, interpolated curve is determined by linear interpolation 
from curve to be generated or processed of the first frame and curve to be generated 
or processed of the third frame to deform the interpolated curve thus obtained into 
curve passing through correspondence points detected at the correspondence point 
detecting section. 
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Operations at the respective steps of FIG. 6, i.e., processing that the arithmetic 
processing unit 2 carries out will now be described in detail. 

At step S102, there is carried out an operation to pursue position of reference 
correspondence points in picture images successive on the time axis. As a pursuit 
method of picture image, there may be employed any methods of pursuing points on 
the picture image. Here, explanation will be given with reference to FIG. 7 in 
connection with a method of carrying out pursuit by using absolute value sum of 
differences of picture images, which is known as SSD as an example. When picture 
image around reference correspondence point 10 at pursuit source frame F 0 shown in 
FIG. 7 is assumed as pursuit source image 20, pursuit destination images 40 are 
prepared with respect to all coordinates within the region of search range 30 
designated at frame F x which is pursuing at present shown in FIG. 7 to determine by 
calculation analogousness between two images, i.e., pursuit source image 20 1? and 
pursuit destination image 40^ As a method of determining analogousness between 
images, there is used a method of determining absolute value sums of differences of 
respective pixel values within image region with respect to two images 20! and 40 x to 
allow them to be evaluation value to allow coordinates having the lowest evaluation 
value within the region to be pursuit result. 

Explanation will now be described in connection with interpolation calculation 
of shape at step SI 03 will be described below. In FIG. 8, there are shown two curves 
where the numbers of segments are not in correspondence with each other. An 
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example for generating intermediate curve when change from first contour curve A to 
second contour curve B takes place will now be described. In this case, when end 
point, control point and line segment of third-order Bezier curve are assumed to be one 
segment, the first contour curve A consists of five segments Sm and the second 
contour curve B consists of three segments Sm. In addition, when intermediate curve 
in the case where change from first contour curve A to second contour curve B takes 
place as shown in FIG. 9, for example, is generated (indicated as intermediate curve 
group Z c to be generated in FIG. 9), procedure shown in the flowchart of FIG. 10 is 
executed. 

When the intermediate curve is generated, correspondence points indicating 
corresponding positions of the first contour curve A and the second contour curve B 
are designated by operation of keyboard, etc. by user. Thus, at step SI, 
correspondence points are set at the first contour curve A and the second contour 
curve B in accordance with the operation to the effect that the above-described 
correspondence points are designated to develop the first contour curve A and the 
second contour curve B into one straight line as shown in FIG. 11 to determine 
respective entire lengths of the first and second contour curves A and B which are 
closed curve to allow the respective entire lengths of the first and second contour 
curves A and B to be length A and length B, respectively. 

Then, at step S2, sampling intervals of the entire length A of the first contour 
curve A and the entire length B of the second contour curve B are determined. 
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Namely, such an approach is employed to compare magnitudes (sizes) of the length 
A and the length B to divide one longer in length by constant of the sampling interval 
to determine the number of points within the section (step S 1 1 of the flowchart shown 
in FIG. 12). In this case 5 constant of the sampling interval is sampling interval where 
there can be obtained the number of points necessary for processing for re- 
constructing curve which will be described later at processing for generating such 
intermediate shape. Then, such an approach is employed to divide smaller one of 
magnitudes of length A and length B by the number of sampling points determined at 
the step Sll to obtain sampling interval (step S12 of FIG. 12). For this reason, 
sampling operations are carried out by sampling intervals less than the maximum 
sampling interval L MAX with respect to both curves at all times. 

Returning to FIG. 10, at step S3, resampling processing of curve A and curve 
B are carried out at sampling intervals determined at the step S2 to prepare point 
trains. Such state is shown in FIG. 13. FIG. 13 shows the state where resampling 
processing are carried out with respect to the first contour curve A and the second 
contour curve B by using sampling intervals obtained at step Sll. In this case, points 
on segment of one Bezier curve can be represented as a function of time t as described 
below. 

Q (t ) = ((l- t ) 3 )M+3t(l-t) 2 )N+3(t 2 )(l-t)0+(t 3 )P (0 s t s 1) 
In the above formula, M, N, O, P are coordinates indicating points on the two- 
dimensional plane, wherein M, P are end points of line segment, and N, O are control 
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points. Further, in order to represent, by t, points on the Bezier curve constituted by 
connecting plural segments, definition of the Bezier curve is extended to allow the 
integral part of t to be segment No. and to allow the decimal part of t to be time within 
segment. Namely, in the case of t=3.45, this indicates coordinate on the curve of the 
third segment and at time of 0.45. Thus, points on the Bezier curve constituted by 
plural segments are represented by t. 

Moreover, detailed flowchart of resampling processing carried out at the step 
S3 of FIG. 10 is shown in FIG. 14. When resampling processing is carried out, times 
at correspondence points of the first contour curve A and the second contour curve B 
corresponding to correspondence point O shown in FIG. 13 are set as Ta and Tb. At 
the subsequent step S22, coordinates of points on the first contour curve A and the 
second contour curve B at times Ta, Tb are determined. Respective points determined 
here result in sampling points at the first contour curve A and the second contour curve 
B. At the subsequent step S23, sampling interval detennined at the above-described 
step S2 of FIG. 10 is added to times Ta, Tb. At the subsequent step S24, there is 
judged whether or not times Ta, Tb reach terminating ends so that correspondence 
points return to correspondence points at the time of start of resampling processing. 
Further, when it is judged that correspondence point returns to correspondence point 
at the time of start, the processing is completed. On the other hand, when it is judged 
that correspondence point does not return to the correspondence point, the processing 
returns to step S22 to repeat processing of the above-described steps S22 to S24 to set, 
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for a second time, sampling points at the first contour curve A and the second contour 
curve B to generate point train consisting of sampling points. Namely, in the 
resampling processing of the step S3 of FIG. 10, there are generated point train 
corresponding to the first contour curve A and point trains corresponding to the second 
contour curve B which consist of sampling points finer than those at the above- 
described sampling processing. 

Further, returning to FIG. 10, at step S4, such an approach is employed to carry 
out interpolation calculation of coordinates of point trains on two curves prepared at 
the previously mentioned step S3 to thereby determine point train of intermediate 
shape. In more detail, processing shown in flowchart of FIG. 15 is carried out to 
thereby generate intermediate curve indicating intermediate shape when there takes 
place change from the contour shape that the first contour curve A indicates to contour 
shape that the second contour curve B indicates. In accordance with the flowchart of 
FIG. 15, first, at step S31, point train corresponding to the first contour curve A and 
point train corresponding to the second contour curve B which are determined by 
resampling processing of the step S3 of FIG. 10 are inputted from the data memory 4. 
Further, there are designated sampling points at positions corresponding to 
correspondence points on point train corresponding to the first contour curve A and 
point train corresponding to the second contour curve B from the correspondence 
points. At the subsequent step S32, when time at which interpolation of two curves 
is carried out is assumed to be T with respect to point train consisting of sampling 
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points, coordinate of point on the first contour curve A is assumed to be a, and 
coordinate of point on the second contour curve B corresponding to the coordinate a 
is assumed to be b, coordinate C of the point at intermediate curve indicating the 
intermediate shape to be determined is expressed by the following formula (2). 



Namely, by applying the above-mentioned formula (2) to point trains corresponding 
to first and second contour curves, coordinate of sampling points constituting 
intermediate curve is determined from sampling points corresponding to the first 
contour curve A and sampling points corresponding to the second contour curve B 
corresponding to sampling points corresponding to the first contour curve A. At the 
subsequent step S33, processing for designating next sampling points of the sampling 
points where the formula (2) is applied at the above-described step S32 is carried out 
with respect to point train corresponding the first contour curve A and point train 
corresponding to the second contour curve B. At the subsequent step S34, the 
arithmetic processing unit 2 judges whether or not correspondence point returns to the 
correspondence point where generation processing of intermediate curve is started at 
the above-described step S31 as the result of the fact that next sampling point is 
designated at the above-described step S33. Further, when it is judged that 
correspondence point returns to the corresponding point, the processing proceeds to 
step S35. On the other hand, when it is judged that correspondence points does not 
return to the correspondence point, the processing which has been explained at the step 



C = Ta + (l-T)-B 



(2) 
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S32 is carried out with respect to sampling points designated at the step S33. Namely, 
the processing which has been explained at the step S32 is carried out with respect to 
all of sampling points set at the above-described sampling processing of the step S3 
to thereby determine coordinate of points at the intermediate curve in accordance with 
the above-mentioned formula (2). At step S35, by using sampling points determined 
at the above-described step S3 and coordinate C of point at the intermediate curve 
determined at the step S32, the technology shown in the Japanese Patent Application 
Laid Open No. 164436/1998 publication or "An Algorithm for Automatically Fitting 
Digitized Curves (pp. 612-626) in the "Graphics Gems" (Andrew S. Glassner) is 
applied to thereby transform (convert) point train consisting of sampling points into 
Bezier curve. Thus, curve is reconstructed from point train consisting of sampling 
points corresponding to first contour curve A and point train consisting of sampling 
points corresponding to the second contour curve B to generate intermediate curve 
indicating intermediate shape when change from first contour curve A to second 
contour curve B takes place as shown in FIG. 9. 

The curve generating apparatus 1 for carrying out such an intermediate curve 
generation processing is adapted so that when, as shown in FIG. 8, intermediate curve 
when there takes place change from the first contour curve A to the second contour 
curve B is generated, even if the numbers of segments of the first contour curve A and 
the second contour curve B are different from each other, the first and second contour 
curves are caused to undergo sampling to set the same number of sampling points with 
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respect to the first contour curve A and the second contour curve B to re-construct 
curve. Accordingly, it is possible to generate intermediate curve group Zc as shown 
in FIG. 9 with high quality and simply. 

Then, at step SI 04 of FIG. 6, interpolated shape determined at step S 103 is 
deformed in correspondence with reference correspondence points determined at the 
step SI 02. Explanation will be given with reference to FIGS. 16, 17 and 18 with 
respect to such an operation. FIG. 16 is the detail of deformation processing of the 
step S 104, FIG. 17 is a view used for explanation from the start START up to the step 
S42 thereof, and FIG. 18 is a view used for explanation from the step S43 up to the 
end END. 

First, at step S41 of FIG. 16, start point 115 and end point 1 16 are caused to 
correspond to each other with respect to reference correspondence points 110 
indicated at START of FIG. 17 determined by image pursuit processing of the step 
S102 of FIG. 6 and reference correspondence points 100 indicated at START of FIG. 
17 determined at the intermediate shape preparation processing of the step SI 03. It 
is now assumed that coordinates of start points at the step S103 are A (exl, ey 1) and 
B(ex2, ey2) and coordinates of end points at the step S102 are C(sxl, syl) and D(sx2, 
sy2). At the step S42, such a coordinate transformation (conversion) to convert line 
segment AB into line segment CD is determined. Such a transformation can be 
represented by combination of parallel displacement and rotation, and operation of 
scaling. Such a transformation is called affine transformation 120. 
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The transform formular of affine transformation is as follows. 
With respect to the method of transforming point (x, y) into point (x\ y') by 
the affine transformation, the following relational expression is assumed to hold, 
[x'] = sc* [cos0 - sin9][x-exl] + [sxl] 
[y'] = sc* [sinG - cos9][y-eyl] + [syl] 
In the above formula, sc = L1/L2 
LI = ^((sx2-sxl) A 2 + (sy2-syl) A 2) 
L2 = N/((ex2-exl) A 2 + (ey2-eyl) A 2) 
6=62-61 

61 = a tan 2 ((sy2-syl)/Ll, (sx2-sxl)/Ll) 
92 = a tan 2 ((ey2-eyl)/Ll, (ex2-exl)/L2) 

Then, at the step S43, affine transformation determined at the step S42 is 
calculated with respect to point train constituting intermediate shape determined at the 
step SI 02 of FIG. 6. As a result, as shown in FIG. 18, there is prepared shape which 
holds feature of intermediate shape determined at the step SI 03 and passing through 
end point determined at the step 102. Then, at step S44, by using the method of 
transforming (converting) point train determined at step S43 into the previously 
described Bezier curve, transformation into curve is carried out. By carrying out such 
a processing with respect to all sections of shape, it becomes possible to transform 
shape prepared at the intermediate shape preparation processing of the step S 103 into 
shape passing through reference correspondence points prepared at the step S102. 
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As described above, in accordance with this invention, in preparing intermediate 
shape of two curves, it becomes possible to prepare intermediate shape passing 
through designated reference correspondence points. In addition, it is possible to set 
reference correspondence point at precise position by image pursuit processing. As 
a result, it has become possible to prepare, with high precision, intermediate shape 
determined at picture image designated on the time axis. Thus, positional shift 
between shape by linear interpolation and actually determined shape which was the 
problem of the prior art can be solved. 

Another more practical example of this invention will now be described. In 
FIG. 19, there is shown the example in which the step SI 03 is omitted among the 
condition of the first embodiment shown in FIG. 6 and the portion of the step SI 04 is 
realized by another method as step SI 05. 

Edge information of picture image is used as substitutive means for preparing 
intermediate shape to carry out processing for automatically generating shape passing 
through the reference correspondence point and along the edge. In more practical 
sense, as a method of detecting optimum path passing through the designated point, 
there is used the technology of path search method disclosed in the Japanese Patent 
Application No. 204862/1997 by the applicant of this application. 

Explanation will be first given in connection with curve generating means which 
generates curve by the path search method disclosed in the Japanese Patent 
Application No. 204862/1997. 
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FIG. 20 is a block diagram showing curve generating means. The curve 
generating means is composed of respective processing means composed of two point 
selector means 15 1, path search processing optimization means 152, path search 
means 153 and curve approximate means 154 5 and respective data stores (storages) 16, 
17, 23, 18 for point train, picture image, path and curve. 

The two point selector means 151 selects successive two points P l5 P 2 from the 
point train data storage means 16 to output them. The path search processing 
optimization means 152 calculates path search parameter P rl for optimizing the path 
search method from the selected two points P l5 P 2 and picture image data D v from the 
picture data storage means 17 to send it to the path search means 153. 

The path search means 153 determines path search processing method on the 
basis of the path search parameter P rl that the path search processing optimization 
means 152 has calculated to calculate path data connecting the selected two points P l9 
P 2 and passing on the contour within picture image to hold it by the path data storage 
means 23. 

In this case, path data is list of pixel coordinates indicating eight vicinity 
graphics connecting the above-mentioned two points. As indicated by broken lines in 
FIG. 20, the path search means 153 sends control information I CN to the two point 
selector means 151 to select the next two points. Accordingly, path searches are 
carried out with respect to portions between respective points of the point train storage 
means 16 so that all the search data are accumulated (stored) into the path data storage 
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means 23. 

When the two point selector means 151 has selected until the last of the point 
train to fail to output new two points, it sends control information I CN to the curve 
approximate means 1 54 as indicated by broken lines in FIG. 20 to hasten the operation 
of the curve approximate means 154. The curve approximate means 154 carries out 
curve approximation of path stored in the path data storage means 23 to prepare curve 
data indicating contour shape to hold it by the curve data storage means 18. 

Explanation will now be given with reference to FIGS. 2 1 and 20 in connection 
with path search method that the curve generating means executes. 

FIG. 21 is a flowchart showing fundamental processing procedure of the path 
search method according to this invention. First, at step SI, the two point selector 
means 151 selects successive two points from the point train data storage means 16 to 
output them. 

Then, at step S2, optimization of path search parameter is carried out. Path 
search parameter P rl is calculated by the path search processing optimization means 
152 on the basis of the two points P 1? P 2 selected by the two point selector means 151 
and picture data Dv. 

Then, at step S3, path search is carried out. The path search means 153 
determines the path search method on the basis of path search parameter obtained at 
the second step S2 to extract contour connecting the above-mentioned two points in 
the form of 8 vicinity path. 
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Further, at step S4, there is carried out judgment as to whether or not path 
searches between all relaying (repeating) points have been completed. As a result, if 
such path searches have not yet been completed, the processing returns to step SI to 
select two points in which one point is shifted from the two points selected at last time 
to repeat steps S2, S3. On the other hand, when it is judged at step S4 that path 
searches between all relaying points have been completed, the processing proceeds to 
step S5. 

At the step S5, contour curve which approximates 8 vicinity path is generated 
to complete the processing. The curve approximate means 154 carries out curve 
approximation of path data stored in the path data storage means 23 to prepare curve 
data indicating contour shape. 

In this case, the curve approximate means 1 54 can be realized by the technology 
in which the applicant of this application has disclosed "Curve Generating Apparatus 
and Curve Generating Method" in the Japanese Patent Application No. 303687/1996 
publication. This technology is a method of obtaining curve which approximates 
coordinate list and is one of the so-called curve approximation methods. Namely, this 
method is method of approximating order of input data while dividing it within the 
sufficiently approximatable range so that this method is robust to noise and has good 
efficiency as compared to the conventional curve approximations. 

However, the above-mentioned curve approximation means 154 may be 
methods of curve-approximating a series of coordinate lists, and may be realized by 
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a method of generating B-spline curve generally used, etc. 

It is to be noted that general curve approximation/interpolation methods are 
described in detail in articles such as "Numerical Calculation Handbook" (Yutaka 
Ohno, Kazuo Isoda, Ohm Company, 1990), "Shape Processing Engineering [I] [II] by 
Computer Display" (Fujio Yamaguchi, Nikkan Kogyo Shinbun Sha, 1982, "Computer 
Graphics Principles and Practice" (Foly, vanDam, Feiner, and Hughes, Addison 
Wesley Company, 1990), etc. 

The path search processing optimization means 152 and the path search means 
153 will now be described in more detail. 

FIG. 22 is a block diagram in which the portion for carrying out path search is 
extracted of the curve generating means. The feature of the curve generating means 
for carrying out this path search resides in that path search processing optimization 
means 152 is provided to thereby carry out path search processing optimized with 
respect to inputted two points. In this example, the path search processing 
optimization means 152 is adapted to output two kinds of path search parameters P r2 
and P r3 which will be described later. 

FIG. 23 shows the internal configuration of the path search processing 
optimization means 152. The path search processing optimization means 152 is 
composed of passing cost calculation range determining means 156 and passing cost 
calculation optimization means 157. The passing cost calculation range determining 
means 156 inputs selected two points P l5 P 2 to output calculation range parameter P r3 . 
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Further, the passing cost calculation optimization means 157 outputs passing cost 
calculation parameter P r2 from the selected two points P l5 P 2 and picture data D v . In 
addition, these two kinds of path search parameters P r2 and P r3 are inputted to the path 
search means 153 and are used for determining the path search method. 

FIG. 24 is a flowchart showing processing procedure at the path search 
processing optimization means 152. 

First, at step SI 1, optimization of passing cost calculation paraneter is carried 
out. Namely, the passing cost calculation optimization means 157 calculates passing 
cost calculation parameter P r2 from the selected two points (point Pj and point P 2 ) and 
picture image data D v . Then, at step S12, optimization of the passing cost calculation 
range parameter is carried out to complete the processing. Namely, the passing cost 
calculation range determining means 156 inputs selected two points P } , P 2 to calculate 
passing cost calculation range parameter P r3 . 

FIG. 25 is a block diagram showing, in more detail, the portion for carrying out 
processing relating to passing cost calculation parameter P r2 in FIG. 23. 

The passing cost calculation optimization means 157 is caused to be of 
configuration comprising gradient parameter determining means 158. 

Moreover, the path search means 153 is caused to be of configuration 
comprising passing cost calculation means 159, minimum cost path calculating means 
161, and passing cost map storage means 160 for storing data called passage cost map 
which will be described later. 
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The gradient parameter deteniiining means 158 inputs picture data D v and two 
points P l9 P 2 to calculate gradient parameter P r4 (a sort of passing cost calculation 
parameters) optimum for extracting contour between the two points to output it to the 
passing cost calculating means 159 of the path search means 153. 

The passing cost calculating means 159 inputs picture data D V3 the point P 1? the 
gradient parameter P r4 , and calculation range parameter P r3 that the passing cost 
calculation range determining means 156 shown in FIG. 23 outputs to calculate 
passing cost map to output it to the passing cost map storage means 160. 

The minimum cost path calculating means 161 inputs passing cost map from the 
passing cost map storage means and the point P 2 to calculate path data D p connecting 
the points Pj and P 2 so as to trace the contour therebetween. 

FIG. 26 is a flowchart showing processing procedure at the path search means 

153. 

First, at step S2 1 , the passing cost calculating means 159 calculates passing cost 
map. Then, at step S22, the minimum cost path calculating means 161 calculates the 
minimum cost path. The detail of this operation will be described later. 

FIG. 27 shows, in more detail, the configuration of the portion for carrying out 
processing relating to the above-described passing cost calculation parameter. 

The passing cost calculation parameter P r2 is two kinds of parameters necessary 
in calculating gradient of picture image. The first parameter is normal vector V 2 
indicating direction perpendicular to advancing direction of contour, and the second 
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parameter is pixel value color change vector V 0 indicating change of pixel value when 
observed in a direction traversing the contour. 

Moreover, the path search means 153 carries out a processing to extract, by 
utilizing the fact that object contour has high gradient strength, such a path on the 
picture image to pass as close as possible the portion having higher gradient strength. 

In order to carry out processing as described above, the pass cost calculating 
means 159 is composed of gradient calculating means 164, passing cost map 
calculating means 165 and gradient data storage means 21. 

The gradient calculating means 164 calculates gradient of picture image. The 
passing cost map calculating means 165 calculates passing cost map such that when 
path is passed through the portion where gradient strength is higher, cost is reduced 
to more degree on the basis of the gradient thus obtained. Position of the point P x is 
inputted to the passing cost map calculating means 165 and starting point of the 
passing cost map is set at the point P 2 . 

In this example, the gradient calculating means 164 selectively detects only 
gradient component of necessary contour by using the technology in which the 
applicent of this application has disclosed in "Edge Detection Method and Edge 
Detection Apparatus" in the Japanese Patent Application No. 134420/1995 
publication. By this gradient calculating means 164, gradient is calculated so as to 
selectively extract gradient of contour extending in a manner connecting inputted two 
points. 
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Explanation will now be given in detail in connection with calculation process 
of the above-described two kinds of path search parameters. 

First, calculation of passing cost calculation parameter P r2 and the process 
where such parameter is utilized will be described. FIG. 28 is a flowchart showing 
processing procedure at the gradient parameter detennining means 158. At step S3 1, 
normal vector calculating means 162 calculates normal vector V x perpendicular to the 
direction connecting inputted points P } and P 2 to output it. At step S32, pixel value 
change vector calculating means 163 calculates change of pixel value of picture image 
when observed in a direction perpendicular to the direction connecting two points, e.g. , 
color change direction in the neighbouring region of the points P L and P 2 on the picture 
image to output it as color change vector V 0 . 

In this example, the passing cost map is information in which,in each pixel, 
direction admitted to that pixel position at minimum passing cost is recorded. 

FIG. 29 shows an example of simple passing cost map. 

The passing cost map is prepared with respect to one range within picture image 
designated by the passing cost calculation range parameter P r3 . In this example, there 
is shown an example where such map is prepared in a certain small quadrilateral 
region. 

In FIG. 29A, directions admitted to that pixel position at the minimum passing 
cost in the respective pixels are indicated by arrow. This direction is calculated so that 
only one direction is determined in each pixel. Value labeled at each pixel is 
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accumulated value of passing cost when admitted in accordance with the minimum 
passing cost direction. This accumulated value is determined by allowing point 
POINT 1 serving as starting point of path search as accumulated value 0. 

FIG. 29B shows process where minimum path calcualting means 161 
determines path when passing cost map is obtained as shown in FIG. 29A. 

The minimum cost calculating means 161 inputs passing cost map and the point 
POINT 2 to trace, in an opposite direction, minimum passing cost direction of the 
passing cost map from the point POINT 2 to thereby determine path of the minimum 
cost ranging from the point POINT 1 to the point POINT 2. For example, in the case 
where the point POINT 2 is located at position shown in FIG. 29 (accumulated cost 
is 22), the arrows are traced from the point POINT 2 (written in a direction opposite 
to that in FIG. 29 A) to determine path until the point POINT 1 as indicated by slanting 
lines. 

Explanation will now be given in connection with means for calculating passing 
cost map as described above at a high speed. Calculation of passing cost map is 
generally known as a sort of optimization problem, and a method of calculating such 
a cost map by using dynamic programming, etc. is well known. Also in the previously 
described article "Intelligent Scissors for Image Composition", a method using 
dynamic programming is described. 

FIG. 30 shows an example of description relating to shortest path search 
algorithm using dynamic programming. It is to be noted that, in this description, signs, 
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etc.. are partially changed from the program described in the above-mentioned article 
"Intelligent Scissors for Image Composition". 

Calculation of calculation range parameter and process where such parameter 
is utilized will now be described. 

FIG. 3 1 is a block diagram showing, in more detail, the portion for carrying out 
processing relating to calculation range parameter in FIG. 23. 

The passing cost calculation range determining means 156 is caused to be of 
configuration comprising distance calculating means 167, calculation range width 
calculating means 168, calculation range length calculating means 169, and calculation 
range calculating means 170. 

The distance calculating means 167 inputs points P l5 P 2 to calculate distance L 
between the two points and direction G connecting the two points. The calculation 
range length calculating means 169 inputs the distance L between the two points 
obtained to calculate calculation range length SL. The calculation range width 
calculating means 168 inputs distance L between the two points obtained to calculate 
calculation range width SW. The calculation range calculating means 170 calculates 
calculation range parameter P r3 to determine shape of the calculation range on the basis 
of the above-mentioned calculation range length SL, calculation range width SW and 
direction G connecting two points that the distance calculating means 167 has 
outputted to input it to the passing cost calculation means 159 of the path search 
means 153. The passing cost calculation means 159 calculates passing cost map 
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within the range indicated by the calculation range parameter P r3 . 

FIG. 32 is a flowchart showing processing procedure at the passing cost 
calculation range determining means 156 of FIG. 31. 

First, at step S41, the distance calculating means 167 inputs points P x and P 2 to 
calculate distance L between these two points and direction G connecting two points. 

Then, at step S42, the distance L between the two points is inputted. Thus, 
calculation range width SW is calculated. In this case, the calculation range width 
calculating means 168 has therein rule of magnitude (size) of calculation range width 
S W with respect to magnitude (size) of distance L between two points in the form such 
as calculation formula or reference table, etc. to determine calculation range width S W 
in accordance with that rule. For example, the problem in which even if operator 
attempts to set relaying (repeating) points in detail at short intervals in order to obtain 
a desired contour, searched path is meandered between relaying points so that shape 
differs from desired shape of contour can be solved by employing a rule such that the 
rule for determining calculation range width S W changes to a rule such that according 
as distance between two points becomes smaller, width becomes smaller. 

At this time, calculation range width SW between relaying (repeating) points set 
in detail at short intervals is set so that it becomes very narrow, and path search is 
carried out so that relaying (repeating) points are connected by substantially straight 
line. Thus, contour shape that operator desires is obtained. 

Then, at step S43, the calculation range length calculating means 169 calculates 




37 



calculation range length SL from distance L between two points. In this case, the 
calculation range length SL is determined so that it is larger than the distance L 
between two points by margin length determined in advance. 

Finally, at step S44, the calculation range calculating means 170 calculates 
calculation range parameter P r3 from the calculation range width S W, the calculation 
range length SL and direction G connecting two points. In this case, the calculation 
range is such a rectangle to surround inputted two points P 1? P 2 and inclined in the 
direction connecting two points, and the calculation range parameter P r3 indicates 
positions of four corners of that rectangle. 

The point train generating means 12 will now be described. 

The point train generating means 12 inputs curve C to output point train 
reproducing that shape and reconstructed curve reconstructed by that point train. The 
point train and the reconstructed curve which have been generated are used for curve 
editing work which will be carried out later by the configuration of the first 
embodiment. 

FIG. 33 is a block diagram showing an example of the configuration of the point 
train generating means 12. 

This point train generating means 12 is caused to be of configuration 
comprising initial point train generating means 171, curve reconstructing means 172, 
difference detecting means 173, point train editing means 14, picture image data 
storage means 17, point train data storage means 16, and reconstructed curve data 
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storage means 19. 

The initial point train generating means 171 inputs input curve C and picture 
data D v to generate initial point train. The point train is stored in the point train data 
storage means 16. The curve reconstructing means 172 inputs point train and picture 
image data D v to generate such a curve to pass the point train and to pass on the 
contour of the picture image. 

The difference detecting means 173 inputs input curve C and re-constructed 
curve C to output difference between both shapes as difference dif The point train 
editing means 14 operates point train data in the point train data storage means 16 on 
the basis of difference dif that the difference detecting means 173 has outputted. 

FIG. 34 is a flowchart showing processing procedure in the point train 
generating means 12. 

First, at step S51, the initial point train generating means 171 carries out 
extraction of feature points of input curve C to store the extracted feature points into 
the point train data storage means 16 as initial point train. 

For the above-mentioned feature point extraction on the curve, there may be 
used feature point extraction method shown in "Motion Vector Estimation Method" 
described in the Japanese Patent Application Laid Open No. 335270/1996 publication. 
In accordance with this method, such an approach is employed to calculate very small 
section bending angles at respective points from the initial end to the terminating end 
of curve to allow its sum total to be accumulated bending angle, whereby when the 
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accumulated bending angle is greater than a threshold value determined in advance, 
the point which takes the maximum very small section bending angle is caused to be 
feature point. Further, similar feature point extraction processing is recursively carried 
out with respect to curve divided by feature points. 

Then, at step S52, the curve re-constructing means 172 calculates re-constructed 
curve by processing similar to the processing in the previously described path search 
apparatus according to this invention from the point train. 

Then, at step S 5 3, the difference detecting means 173 calculates difference dif 
between input curve and re-constructed curve. In accordance with such difference 
detection processing, such an approach is employed to take difference detection points 
at suitable intervals on the re-constructed curve to measure distance up to input curve 
from that difference detection point to allow it to be difference value dif. 

Then, at step S54, the difference detecting means 173 examines whether or not 
there exist difference value dif greater than threshold value set in advance. If there 
does not exist difference value greater than the threshold value, the processing of the 
point train generating means 12 is completed. Thus, data in the point train data storage 
means 16 and re-constructed curve data storage means 19 at that time point serve as 
output. On the other hand, if there exists difference value greater than the threshold 
value, the processing proceeds to step S55. 

At the step S55, the point train editing means 14 edits point train data on the 
basis of difference dif thus obtained. At respective difference detection points in 
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which it is judged that the difference value dif is greater than threshold value, new 
relaying (repeating) point is placed at point on the input curve nearest from that 
difference detection point to store point train data into the point train data storage 
means 16. 

After completion of the step S55, the processing returns to the step S52. Thus, 
at step S54, processing is repeated until difference between input curve C and re- 
constructed curve C becomes smaller than threshold value. 

While, in the above-mentioned processing, the number of points constituting 
initial point train is set to lesser value at first to suitably supplement initial points at the 
shortage portions, there may be employed a method in which the number of points 
constituting initial point train is set to larger value to delete unnecessary points. In that 
case, such an approach may be employed to generate initial point train to such a degree 
that input curve can be sufficiently polygon-approximated by the initial point train 
generating means 171 to select and delete relaying (repeating) points in which 
difference value therearound is sufficiently smaller than threshold value by the 
difference detecting means 173 after re-construction of the curve. 

It is to be noted that while, in the processing procedure shown in FIG. 34, since 
re-construction can be completed by lesser number of points, there results the merit 
that curve editing thereafter is easy, processing procedure which will be described later 
has the merit that because relaying (repeating) point interval for carrying out path 
search can be small, processing time can be shortened. 
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Moreover, point train obtained by the above-mentioned method is transformed 
(converted) into curve by using curve generating means. Namely, by using edge 
information of picture image in preparing intermediate shape passing through portions 
between respective correspondence points by using the above-mentioned method, it 
has become possible to prepare contour having high precision along the shape 
represented in the picture image. 

Further, in the method shown in FIG. 19, such an approach may be employed 
to determine edge of picture image by using another algorithm at step S105 to 
determine intermediate shape along the edge of that picture image. 

Namely, at step S102, a set of correspondence points determined by picture 
image pursuit and the picture image are used to determine color projection axis which 
maximizes contrast around the edge to form edge of picture image by using such color 
projection axis to determine intermediate shape along this edge. 

The configuration of optimum color projection axis calculating unit 32 which 
inputs the set of correspondence points as representative points Ql and Q2 to 
determine the color projection axis P jo is shown in FIG. 35. 

The optimum color projection axis calculating unit 32 is composed of a 
representative point vicinity (neighbouring) color projection axis calculating section 
41a and a representative point vicinity (neighbouring) optimum color projection axis 
calculating section 41b for calculating vicinity (neighbouring) optimum color 
projection axes and P j2 within the region in the vicinity of representative points Ql 
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and Q2, and a color projection axis synthesis section 42 for synthesizing color 
projection axis P j0 by using the color projection axis P }} and the color projection axis 
P j2 from the above-mentioned two calculating sections 41a, 41b. 

At the representative point vicinity (neighbouring) optimum color projection 
axis calculating sections 41a and 41b, there are determined vicinity (neighbouring) 
optimum color projection axes and P j2 which maximize contrast around the edge 
included within neighboring regions of representative points Ql and Q2. It is now 
assumed that, as shown in FIG. 36, representative points Ql and Q2 are determined 
within the regions color-separated by black B, white W, green G, yellow Y 5 red R and 
blue B. At this time, the representative point vicinity (neighbouring) color projection 
axis calculating section 41a determines, by calculation, color projection axis P^ in 
which the contrast is maximum of edge included within the neighboring region Al of 
representative point Ql. In addition, the representative point vicinity neighboring 
color projection axis calculating section 41b determines, by calculation, color 
projection axis P j2 in which the contrast is maximum of edge included within the 
vicinity (neighbouring) region A2 of the representative point Q2. 

The color projection axis synthesis section 42 synthesizes vicinity 
(neighbouring) optimum color projection axis P jl9 vicinity (neighbouring) optimum 
color projection axis P j2 respectively independently determined within the vicinity 
(neighbouring) regions of representative points Ql, Q2 to determine one color 
projection axis P jo where the contrast is maximum of edge of updating section. As the 




43 



synthesis method in this case, e.g., average vector of color projection axes may be axis 
after synthesis. 

In addition, the maximum color projection axis calculating section 32 may 
determine optimum color projection axis by the method shown in FIG. 37. 

As shown in FIG. 3 7 A, it is assumed that anchor point ApO and adjacent anchor 
point Apl are placed on edge e of picture image. At this time, since there is the 
possibility that the edge e may deviate to much degree from the line segment 
connecting ApO and Apl, there is carried out search of region where color change is 
the greatest within the range of length H on the line segment in a direction 
perpendicular to ApO Apl with respect to three points placed every equal interval with 
respect to line segment length L on the line segment ApOApl, e.g., interval of L/4. 

When color with respect to a certain point r is assumed to be C(r), it is defined 
that color change vector is dC/dr and its magnitude is |dC/dr|. Thus, r in which this 
value takes the maximum value is determined. 

As the result of search, rO, r 1 , r2 are assumed to be determined as shown in FIG. 
37B. At this time, color change vectors dC/dr are respectively determined at rO, rl, 
r2 to allow its average value to be optimum color projection axis. 

In addition, color change vector may be determined in a manner described 
below. In the case of rO, colors C(sOO), C(s01) are determined with respect to points 
sOO, sOl having rO as center and away from each other by distance width h on the 
straight line perpendicular to the line segment ApOApl. Thus, color change vector is 
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substituted for its difference vector aC=C(s01)-C(s00). In connection with rl, r2, 
color change vectors are respectively similarly determined. 

Ordinary contour closed curve generating unit generates contour closed curve 
on the basis of the optimum color projection axis which has been explained. At the 
step S105, intermediate shape based on this contour closed curve is prepared. 

In addition, when algorithm shown in FIG. 19 is deformed and linear 
interpolation at the step S 102 fails to be carried out because, e.g., movement of picture 
image is suddenly carried out, the processing of the step SI 05 may be carried out. 

As described above, in accordance with this invention, positional shift between 
shape by linear interpolation and shape actually determined, which was problem of the 
prior art, can be solved. 

Industrial Applicability 

As previously mentioned, in accordance with this invention, even in the case 
where contour shape does not linearly move, it is possible to precisely pursue 
correspondence points on the contour. Thus, contour shape having higher precision 
can be prepared. 



